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ABSTRACT 

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a 
glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide 
meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the 
stars are among the brightest in the Kepler field. We present an analysis of three months of nearly 
uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 
46 and 41 oscillation frequencies for the two components respectively, including a clear detection of 
octupole (1=3) modes in both stars. We derive the properties of each star independently using the 
Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational 
constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results 
generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by 
fitting each component individually yield a common age (t = 6.8 ± 0.4 Gyr) and initial composition 
[Z\ = 0.024 ± 0.002, Yi = 0.25 ± 0.01) within the uncertainties, as expected for the components of a 
binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data 
sets that will ultimately become available will allow future studies of differential rotation, convection 
zone depths, and long-term changes due to stellar activity cycles. 

Subject headings: stars: fundamental parameters — stars: individual(HD 186408, HD 186427) — stars: 
interiors — stars: oscillations — stars: solar-type 



1. INTRODUCTION 

Observations of the Sun provide an extraordinarily de- 
tailed snapshot of stellar structure and dynamics for a 
single set of physical properties and at a particular evo- 
lutionary state in the life of a star. To generalize our un- 
derstanding of stellar evolution, and to evaluate the de- 
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gree to which the Sun is typical or peculiar, it is useful to 
examine other classes of stars. iCavrel de Strobell ([1996H 
defined several such classes, including: solar twins with 
fundamental physical properties very similar or identical 
to the Sun; solar analogs which are broadly comparable 
to the recent past and near future of the Sun; and solar- 
like stars including a wider range of F and G dwarfs 
and subgiants. W ell-known solar twins such as 18 Sco 
(Baz ot et al1l201lH provide some context for the Sun ob- 
served as a star, whi le solar analogs like k 1 Cet, (3 Hyi, 
and a Cen A & B (iWalker et all 120071: iBrandao et all 



IMi^lBedding et al.ll2004l: IKieldsen et~aLll2005D help cal- 
ibrate stellar evolution for stars that are younger, older, 
and more or less massive than the Sun. Broader studies 
of solar-like stars probe the full range o f relevant stel- 
lar properties and evolutionary states (e.g. lChaplin et all 
IMllSilva Aeuirre et alJ[201ll) . 

The bright stars 16 Cyg A & B (= HD 186408 & 186427 
= KIC 12069424 & 12069449; V~6) have long been stud- 
ied as solar analogs , with estimated ages n ear 6-8 Gyr 
([Wright et al.1 120041 : IVaTenti fe Fischer! 12005( 1 . Although 
they are members of a hierarchical triple system with 
a red dwarf companion that is 10 mag nitudes fainter 
(|Turner et all 120011: iPatience et all 12002( 1. there are no 
dynamical constraints on the masses because the avail- 
able d ata suggest an orbital period longer than 18,000 
years ([Hauser fe Marcyl 11999( 1. After the discovery of a 
1.5 Jupiter- mass exop lanet in an eccentric 800-day or- 
bit around 16 Cyg B (|Cochran et al.lll997l) . the system 
generated even more interest. Since then, both compo- 
nents have been monitored for magnetic activity, showing 
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Figure 1. Power spectra of 16 Cyg A (left panels) and 16 Cyg B (right panels). Top panels: 20 fiHz boxcar smoothed spectra (grey), 
with best-fitting background components attributed to granulation (dashed lines), stellar activity and/or larger scales of granulation (dot- 
dashed lines) and shot noise (dotted lines), with the sum of the background components plotted as solid black lines. Bottom panels: 
Background-subtracted power spectra over the ranges in frequency where high-order p modes are observed. 



long-term variations around a mean chromospheric ac- 
tivity level well below that of the Sun at solar minimum 
(J. Hall, private communication). So far there have been 
no direct measurements of rotation, but gyrochronology 
suggests that the r otation periods should be near 30 days 
(Skumanich 1972). As two of the brightest stars in the 
Kepler held of view, 16 Cyg A & B can now be subjected 
to a long-term study that promises to yield more detailed 
information than is currently available for any star but 
the Sun. 

In this Letter, we perform an asteroseismic analysis of 
the first three mont hs of data on 16 Cyg A & B from 
the Kepler mission (jKoch et al.1 120101 ) . Using the un- 
precedented observations, we model the two components 
independently and determine an identical age and initial 
composition, as expected for the members of a binary 
system. In $2] we describe the data analysis methods, 
and in £j3]we present the asteroseismic modeling includ- 
ing an evaluation of both the statistical and systematic 
uncertainties. We conclude in $4] with a discussion of the 
results and the potential for future studies utilizing the 
longer data sets that will soon become available. 

2. DATA ANALYSIS 

The first full q uarter of short-cade nce observations 
(58.85 s sampling; IGilliland et all[20l0l) of 16 Cyg A & 
B were obtained by Kepler between September and De- 



cember 2010 (Q7). Both stars are significantly brighter 
than the photometric saturation limit. Saturated flux is 
conserved on Kepler, so no photometric precision is lost 
for saturated targets as long as the saturated pixels are 
included in the pixel aperture. Standard Kepler pixel 
apertures were not designed for such bright, saturated 
targets and in the case of 16 Cyg contain a prohibitively 
large number of unneeded pixels. Custom masks were 
therefore defined from Q3 full frame images to capture 
all of the flux using fewer pixels. The ra w photometric 
light- curves extracted from these masks (|Jenkins et al.l 
120101 ) were then prepare d for asterose i smic analysis in 
the manner described by iGarcia et al.l (|2Qlll ). Figure Q] 
shows the power spectra of both stars (16 Cyg A in the 
left panels, and 16 Cyg B in the right panels). 

The top panels include boxcar smoothed power spectra 
(in grey) over an extended range in frequency, showing 
not only the Gaussian-like power excess due to solar- 
like oscillations, but also contributions to the back- 
ground power-spectral density attributable to granula- 
tion (dashed lines), stellar activity and/or larger scales 
of granulation (dot-dashed lines) and shot noise (dotted 
lines) . The backgrounds were fit with a three-co mponent 
mode l, comprising two Harvey-like power laws (|Harvevl 
1985) to represent granulation and activity, and a flat 
component to represent the contribution of shot noise. 
The best-fitting sum of background components is shown 
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Table 1 

Observed oscillation frequencies for 16 Cyg A 



16 Cyg A 



16 Cyg B 



1 = (pHz) 



= 1 (pHz) 



! = 2 (/iHz) 



1 = 3 (mHz) 



= fuHz) 



= 1 (/iHz) 



= 2 (/iHz) 



= 3 fugz) 



13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



1598. 
1700. 
1802. 
1904. 
2007. 
2110. 
2214. 
2317. 
2420. 
2524. 
2629. 
2736. 
2838. 



51 ± 0.27 
43 ± 0.34 
15 ± 0.17 
62 ± 0.15 
45 ± 0.13 
94 ± 0.11 
33 ± 0.17 
18 ± 0.17 
75 ± 0.30 
94 ± 0.39 
36 ± 0.36 
22 ± 1.45 
68 ± 0.38 



1644 
1746 
1849. 
1951 
2055. 
2158. 
2262. 
2366. 
2470 
2575. 
2678 
2783. 
2889. 



24 ±0.33 
93 ± 0.24 
11 ± 0.13 
98 ± 0.16 
41 ± 0.16 
89 ± 0.12 
32 ± 0.16 
15 ± 0.16 
23 ± 0.25 
97 ± 0.31 
47 ± 0.47 
71 ± 1.22 
61 ± 0.38 



1591 
1693 
1795 
1898 
2001 
2105 
2208 
2312 
2416 
2520 
2624 
2730 



21 ± 0.86 
73 ± 0.46 
87 ± 0.40 
.08 ± 0.27 
± 0.17 
.60 ± 0.15 
.90 ± 0.19 
.49 ± 0.29 
.24 ± 0.33 
91 ± 0.81 
.05 ± 0.51 
06 ± 1.03 



1736. 
1839. 
1944. 
2045. 
2150. 
2253. 
2356. 
2461. 



03 ± 1.84 
07 ± 1.64 
07 ± 1.57 
09 ± 0.80 
15 ± 0.19 
41 ± 0.35 
92 ± 0.46 
26 ± 1.04 



1928 
2044 
2159 
2276 
2392 
2509 
2626 
2743 
2860 
2978 
3096 
3215 



1 ± 0.28 
.21 ± 0.15 
36 ± 0.16 
.03 ± 0.12 
.87 ± 0.14 
.75 ± 0.13 
43 ± 0.11 
15 ± 0.25 
63 ± 0.26 
95 ± 0.40 
.00 ± 0.54 
.94 ±0.91 



1982 
2098 
2214 
2330 
2448 
2565 
2682 
2799 
2917 



66 ± 0.16 
20 ± 0.17 
00 ± 0.18 
88 ± 0.16 
17 ± 0.11 
35 ± 0.10 
38 ± 0.14 

67 ± 0.22 
75 ± 0.22 



1920. 
2036. 
2152. 
2269. 
2386. 
2503. 
2619. 
2737. 
2854. 
2972. 
3089. 



99 ± 0.24 
59 ± 0.20 
91 ± 0.19 
07 ± 0.21 
30 ± 0.17 
56 ± 0.13 
99 ± 0.23 
44 ± 0.31 
52 ± 0.39 
73 ± 0.70 
46 ± 0.87 



2202.75 ± 0.65 
2317.08 ± 0.44 
2436.78 ± 0.33 
2553.00 ± 0.23 
2672.34 ± 0.28 
2788.74 ± 1.40 
2906.96 ± 0.93 



3152.45 ± 0.61 
3274.63 ± 0.55 



i Radial order n from the optimal AMP models. 



as a solid black line in each panel. We note that the best- 
fitting timescales (T gra n,A = 257 ± 6 sec, T gra n,B = 241 ± 
8 sec) and peak powers (P gra n,A = 3.01±0.08 ppm 2 /^Hz, 
-Fgran.B = 2.41 ±0.07 ppm 2 /^Hz) of the granulation com- 
ponents are both slightly greater than the solar values es- 
timated from Sun-as-a-star observations, and fol l ow th e 
scaling relations derived by Kicldscn & Bedding ( 20111) . 

The bottom panels show very clear patterns of peaks 
due to solar-like oscillations of high radial order, n. The 
quality of the oscillation spectra are exquisite, with each 
star showing more than fifteen radial overtones including 
many octupole (^=3) modes. The maximum peak height- 
to-background ratios are comparable to those observed 
in photometric Sun-as-a-star data — the shot noise level 
is so low in these Kepler data that the intrinsic stellar 
granulation noise dominates the background across the 
frequency ranges where the modes are observed. 

Ten teams provided estimates of the frequencies of the 
observed modes, applying peak-baq qinq techniques de 



velop ed for application to CoRoT (lAppourchaux et al 



2008) 



2011 



and Kepler data (e.g., see ICampante et~ 
IMathur et al.l 1201 ID . These techniques varied 
in the details of the optimizations performed — which 
included class i cal m aximum-likelihood estimation (e.g. 



iFletcher et al 



120111) 



and Markov Chain Monte Carlo 
methods (e.g. lHandbcr g fc Camp antc 201l|) — and in the 
number of free parameters and assumptions made for fits 
of Lorentzian-like models to mode peaks in the power 
spectra. 

The results of the ten teams were analyzed to produce 
final frequency sets for each star. First, we sought to 
identify objectively those modes for which a robust, well- 
determined frequency could be estimated. This involved 
two types of checks. In one, we identified a list of modes 
with good agreement between the best-fitting frequencies 
from various teams. This was achieve d using a modi- 
fied ve rsion of the procedu r e outl ined in ICampante et al.l 
(|201in and IMathur et al.l (|2011D . A so-called minimal 
frequency set of modes was produced, for which a major- 
ity of the teams' estimates were retained after applying 
Peirce's criterion ([Peir cdl 1 8 5 2t I Gouldll 1 8 5 51 ) for outlier re- 
jection. A second set of checks involved visual inspection 
of the frequency- power spectra (and echelle diagrams of 
those spect ra), combined with o bjective false alarm prob- 
ability (e.g. lChaplin et al1l2002D and likelihood ratio tests 



(e.g. I Appourchauxl 120 1 ID . 

With a list of robust modes in hand, one of the teams 
was then selected to re-fit these modes in both stars us- 
ing a single Lorentzian profile per mode (i.e. no rotational 
splitting and the inclination angle fixed at 0°). This team 
was chosen as the one whose initial best-fitting frequen- 
cies showed the closest match to the frequencies of the 
minimal-set. Use of frequencies from one of the teams, 
as opposed to some average over all teams, meant that 
the modeling could rely on an easily reproducible set of 
input frequencies (see Table [IJ. 

3. ASTEROSEISMIC MODELING 

The set of oscillation modes from the peak-bagging 
analysis described in [21 included a total of 46 and 41 
individual frequencies for 16 Cyg A & B, respectively. 
As inputs for the stellar modeling, we supplemented 
these asteroseismic constraints with the spectroscopic 
prope rties of each component derived by iRamirez et al.l 
(2009) [T off A = 5825 ± 50 K, log 5A = 4.33 ± 0.07, 
[Fe/TfjA =' 0.096 ± 0.026; T eS B = 5750 ± 50 K, 
log(? B = 4.34 ±0.07, [Fe/H] B = 0.052 ± 0.021]. Us- 
ing t hese T fi ff valu es to obtain bolometric corrections 
from iFlowerj iHHHi: and adopting M hol = 4.73 ± 0.03 
from lTorresllhoiOD. we combined the extinction estimates 
from lAmmons et al.l (120061) wi th the updated Hipparcos 
parallaxes (|van Leeuwenll2007l ) to derive luminosity con- 
straints: L A = 1.56 ± 0.05 L Q , L B = 1.27 ± 0.04 L Q . 

We calculated separate values of x 2 f° r the astero- 
seismic and spectroscopic constraints and attempted to 
minimize them simulta neously using the A steroseismic 
Modeling Portal (AMP; lMetcalfe et all l2009l). This auto- 
mated method uses a parallel genetic algorithm to search 
a broad range of stellar parameters and objectively de- 
termines the globally optimal model for a given set of 
observations. Although 16 Cyg A & B are members of 
a binary system and presumably formed simultaneously 
from the same material, we fit each set of constraints 
independently and did not force the models to have a 
common age or initial composition. The oscillation fre- 
quencies of the optimal models for each star are plotted 
as red symbols in Figure [H where the observed modes 
are shown as black points with horizontal error bars. In 
both cases, the asteroseismic x 2 is less than 10 and the 
spectroscopic \ 2 is l ess than 1, so the models represent 
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Figure 2. Bchclle diagrams of 16 Cyg A (left) and 16 Cyg B (right), showing the observed frequencies as black points with horizontal 
error bars. The frequencies of the optimal models from AMP are shown using different red symbols to indicate modes with radial (1=0, 
circles), dipole (1=1, triangles), quadrupole (1=2, squares) and octupole (1=3, diamonds) geometry. A greyscale map showing a Gaussian 
smoothed power spectrum (FWHM ~ 2 /iHz) is included in the background for reference. 



a reasonably good match to both sets of observational 
constraints. In Table [5] we list the optimal values for the 
radius (R) and for the adjustable model parameters, in- 
cluding: the mass (M), age (t), initial metallicity (Z\) 
and helium mass fraction (Y\), and the mixing- length 
parameter (a), along with the asteroseismic % 2 . The 
statistical uncertainties on each parameter (<7 s t a t) were 
determined using Singular Value Decomposition. 

To evaluate the possible sources of systematic uncer- 
tainty from the ingredients and assumptions in our mod- 
els, six teams were given the results from AMP and asked 
to reproduce the fit using the same set of observational 
constraints with their own stellar evolution codes and fit- 
ting methods. The physical ingredients adopted by each 
team differed slightly from those employed by AMPr^l. 
allowing us to explore the degree of model-dependence 
in our results. The optimal parameter values from each 
team are listed in Table [2j where we combine all of them 
into an adopted value (bold row) representing the aver- 
age of the individual estimates weighted by 1/x 2 - The 
systematic uncertainty (cr sys ) on each parameter reflects 
the variance of the results, again weighted by 1/x 2 - 

As expected, there are slight differences between the 
optimal parameter values determined by each team. 
Since we effectively used AMP to solve the global op- 



15 AMP uses the OPAL 2005 equation of state and the most 
recent OPAL opacities supplemented by Alexander & F erguson] 
1119941) opacities at low tempe rature, nuclear reaction rates from 
Bahcall & Pinsonncault (19921), and includes the eff e cts of helium 
diffusion and settling following Michaud & Proffitt (1993). Con- 
vection is treated with standard mixing-length theory without over- 
shooting (Bohm-Vitcnsc 1958). ANKI so lves the Saha e q uatio n 
itself, uses low temperature opacities from IFereuson et alj (2005), 
nuclear reaction rates primarily from Caughlan & Fowler (1988), 
and in cludes a full treatment of di ffusion followin g ITh oul et al.l 
(IT99D . ASTEC1 and AST EC2 use [F erguson et al. (2005) opaci- 
ties, NACRE reaction rates (Angulo et al. 1999), and neglect dif- 
fusion. CESA M uses NACRE rea ction rates and treats convec- 
tion following ICa nuto et al. (1993). Gen eva uses NACRE reac - 
tion rates, and treats diffusion f ollowing Paquettc e t al.l (119861 ). 
YREC uses Ferguson et al. ( 2005) opac ities, nuclear reaction rates 
primarily from [A dclbcrgcr et al. (1998j), treats diffusion following 
IThoul et al. ( 1994) and includes convective overshooting. All mod- 
els include the empirical correction for surface effects proposed by 



els include the emp: 
IKieldsen et~aTI (2008). 



timization problem, these differences reflect subtle shifts 
in the locally optimal solution due to the physical ingredi- 
ents included in each stellar evolution code. However, the 
results from different teams also include small offsets due 
to incomplete optimization — refined sampling of each ad- 
justable parameter will always improve the fit, and there 
was no uniform criterion for when to stop fine tuning. To 
minimize the influence of this technique error on the final 
results, we weight the average parameter values and un- 
certainties using 1/x 2 from each result as a proxy for the 
overall quality of the fit. This ensures that the variance 
reflects the actual systematic differences between model 
physics rather than the effort expended by each team in 
trying to match the observations. As with AMP, most 
of the teams did not force any of the model parameters 
of 16 Cyg A & B to share a common value. The ex- 
ceptions were the Geneva code (which forced a common 
age, initial composition and mixing- length), and YREC 
(which forced the model for B to have the same age as 
the optimal model for A). Excluding these models from 
the average does not significantly change the values of 
the adopted parameters listed in Table dJ 

4. RESULTS & DISCUSSION 

We have performed an analysis of the solar analogs 
16 Cyg A & B using three months of observations from 
the Kepler space telescope, yielding the highest quality 
asteroseismic data sets for any star but the Sun (see 
Figured]). We identify a total of 46 and 41 oscillation 
frequencies in the two components respectively, includ- 
ing a clear detection of octupole (1=3) modes in both 
stars. These modes are difficult to detect in photomet- 
ric data because the bright and dark patches associated 
with higher degree modes are normally expected to can- 
cel in disk-integrated measurements. The unambiguous 
detection of such modes from the Kepler light-curves of 
16 Cyg A & B is a testament to the exceptional quality 
of the data. 

We derived the properties of each star independently 
by fitting stellar models to the oscillation frequencies 
(see Table [T]) and other observational constraints (see 
Sj3]) simultaneously. The initial results from fitting each 
star individually using the Asteroseismic Modeling Por- 
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Table 2 

Stellar model-fitting results for f6 Cyg A & B. 



16 Cyg A 16 Cyg B 





R/R @ 


M/Mq 


t(Gyr) 


Zi 


Y 


a 


x 2 


R/R e 


M/Mq 


t(Gyr) 


Zi 


Yi 


a 


x 2 


AMP 


1.236 


1.10 


6.5 


0.022 


0.25 


2.06 


5.47 


1.123 


1.06 


5.8 


0.020 


0.25 


2.05 


9.80 




0.016 


0.01 


0.2 


0.002 


0.01 


0.03 




0.020 


0.01 


0.1 


0.001 


0.01 


0.03 




ANKi 


1.260 


1.14 


6.4 


0.024 


0.26 


1.94 


21.41 


1.138 


1.08 


6.4 


0.022 


0.26 


1.94 


23.29 


ASTEC1 . . 


1.237 


1.10 


7.5 


0.023 


0.25 


2.00 


5.70 


1.121 


1.05 


7.3 


0.021 


0.25 


2.00 


7.97 


ASTEC2 . . 


1.235 


1.10 


6.8 


0.022 


0.25 


2.00 


7.70 


1.134 


1.09 


6.3 


0.025 


0.25 


2.00 


8.47 


CESAM. . . 


1.253 


1.14 


7.0 


0.027 


0.24 a 


0.72 b 


3.53 


1.136 


1.09 


6.9 


0.025 


0.24 a 


0.73° 


4.78 


Geneva. . . . 


1.236 


1.10 


6.7° 


0.024° 


0.26 c 


1.80 c 


10.82 


1.122 


1.06 


6.7° 


0.024° 


0.26 c 


1.80 c 


10.98 


YREC .... 


1.244 


1.11 


6.9 


0.026 


0.26 


2.08 


5.68 


1.121 


1.05 


6.9 d 


0.022 


0.26 


1.84 


3.17 


adopted 


1.243 


1.11 


6.9 


0.024 


0.25 


2.00 




1.127 


1.07 


6.7 


0.023 


0.25 


1.92 






0.008 


0.02 


0.3 


0.002 


0.01 


0.08 




0.007 


0.02 


0.4 


0.002 


0.01 


0.09 





a Values of Yi < 0.24 excluded from search. 

b Value of a from the Canuto ct al. (1996) treatment of convection, excluded from average. 
c Age, composition, and mixing-length constrained to be identical in both components. 
d Age of 16 Cyg B constrained to be identical to the value found for 16 Cyg A. 



tal (AMP; see Figure ^ yield the same initial compo- 
sition within the statistical uncertainties, and a similar 
age for the two components. Further analysis using sev- 
eral stellar evolution codes employing a variety of input 
physics allows us to quantify the model-dependence of 
our results, and to adopt reliable values and uncertain- 
ties from the ensemble. The adopted stellar properties of 
16 Cyg A & B (see Table[2|) reinforce the conclusion that 
the two stars share a common age (t = 6.8±0.4 Gyr) and 
initial composition (Zi = 0.024±0.002, Y ; = 0.25±0.01), 
as expected for a binary system but without imposing 
this as a constraint for the modeling. This fundamental 
result bolsters our confidence in the reliability of astero- 
scismic inferences of stellar structure and evolution. 

The relative size of the statistical and systematic un- 
certainties provides an important benchmark for what we 
can expect from asteroseismology with the Kepler mis- 
sion. The statistical uncertainties on the stellar radii 
from AMP were derived from the distribution of radii in 
an ensemble of models that differ from the optimal model 
by ±1(7 for each adjustable parameter. Such estimates 
implicitly include the influence of parameter correlations, 
and are consequently much larger than the systematic 
variation in optimal radii from different stellar evolution 
codes. This is not the case for adjustable model param- 
eters like the mass and age, where a sys can be 2-4 times 
larger than cr s t a t, while the two are roughly comparable 
for the initial metallicity and helium mass fraction. The 
mixing-length parameter is a special case, because the 
range of estimates from different stellar evolution codes 
reflect variations in the solar-calibrated values that arise 
from differences in the input physics and in the specific 
formulation of mixing-length theory that is implemented 
in each code. Thus, the systematic uncertainties on a 
are likely to be overestimated, and small variations in 
the optimal value of this parameter from different codes 
should not be overinterpreted. 

These extraordinary results were possible using just 
the first three months of short-cadence observations (Q7) 
from Kepler. Nine months of data will soon be avail- 
able (Q7-8-9), and the stars continue to be on the short- 
cadence target list — at least through Q12, and hopefully 
for the remainder of the mission. These longer data 
sets will gradually yield higher frequency precision and 



improve the signal-to-noise ratio in the power spectra, 
enabling further characterization of the stars from the 
current frequency sets and facilitating detection of addi- 
tional oscillation modes at higher and lower frequencies. 

From 6-9 months of data we may begin to resolve ro- 
tational splitting of the non-radial oscillation modes into 
multiple components with different azimuthal order, m. 
The variation of this splitting as a function of the ra- 
dial order n can probe radial differential rotation, while 
the differences between non-radial modes with different 
spherical degree I can reveal latitudinal variations. Such 
measurements of rotation may help to constrain possi- 
ble scenarios to explain the different Li abundances o f 
the two stars (jSchuler et all 120111 : iRamfrez et all 120 111 ). 
With 12-18 months of data, the frequency precision may 
be sufficient to resolve oscillatory signals in the devi- 
ations from uniform frequency spacing (so-called sec- 
ond differences, 82V) which reflect the acoustic depths 
of sharp transitions in the stellar structure, such as the 
helium ionization region and the base of the surface con- 
vection zone. Even longer data sets will allow us to probe 
the influence of stellar activity cycles, which lead to small 
anti-correlated changes in the frequencies and amplitudes 
of the oscillation modes. By the end of the baseline Ke- 
pler mission, and hopefully through an extended mission, 
these two bright solar analogs promise to yield the clear- 
est picture yet of the future of our own Sun. 
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